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Abstract 

Recent work on the connection between in-medium subthreshold K~N 
amplitudes and kaonic atom potentials is updated by using a next to leading 
order chirally motivated coupled channel separable interaction model that 
reproduces KN observables at low energies, including the very recent SID- 
DHARTA results for the atomic i^~-hydrogen Is level shift and width. The 
corresponding ^"-nucleus potential is evaluated self-consistently within a 
single-nucleon approach and is critically reviewed with respect to empirical 
features of phenomenological optical potentials. The need to supplement 
the single-nucleon based approach with multi-nucleon interactions is demon- 
strated by showing that additional empirical absorptive and dispersive terms, 
beyond the reach of chirally motivated .^"-nucleus potentials, are required in 
order to achieve good agreement with the bulk of the data on kaonic atoms. 
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1. Introduction 

The bulk of the data on strong interaction effects in kaonic atoms is due 
to experiments of over three decades ago. With the exception of the very 
light atoms of K~H and K~ A He the rest of the data could be described rather 
well with the help of -^"-nucleus optical potentials [l|,0]. Recent experiments 
on K~H and K~ 4 He with much reduced background removed the 'puzzles' 
with these two atoms j^l, 0, 0, @] . However, the depth of the attractive K~- 
nucleus real potential, which in phenomenological analyses came out in the 
range of 150-200 MeV j?J , presented a theoretical challenge in as much as in- 
medium chiral threshold K~ N scattering amplitude input led to a lower value 
of order 120 MeV jsf, or even to a considerably lower value of 40-50 MeV jof. 
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This outstanding discrepancy is of current interest, since it is relevant to the 
role of K~ mesons in multistrange self-bound matter and in compact stars 
[10 . [Tlj] . The problem has been largely resolved very recently (l^, 13] noting 



that in- medium chiral subthreshold K~N scattering amplitudes provide the 
relevant input and thereby demonstrating the need to supplement the model 
by multi-nucleon terms, as discussed in the present work. This leads to deep 
real potentials in agreement with the purely phenomenological analyses. 



The present paper is an update of Refs. [12j, [13J, based on a recent 
in-medium coupled channel chirally motivated separable interaction model 
which produces good fits to all the low energy antikaon-nucleon data, in- 
cluding the latest fT~H atom results from the SIDDHARTA experiment [if. 
Section [2] outlines the self-consistent handling of subthreshold K~N ampli- 
tudes while section [3] deals with the resulting ^"-nucleus amplitudes. Section 
H] reports on global optical model fits to kaonic atom data and the last section 
provides summary and conclusions. 

2. K~N scattering amplitudes 

The potential experienced by a K~ meson of energy E l ^ h = uk interacting 
with a nucleus of density p is given in the single-nucleon approximation by 

V K - (u K ; p) = (1 + — ) F K ~ N (p, y/s; p) p, (1) 

u K m N 

where Fk-n{p, p) is the in- medium K~N scattering amplitude, reducing 
in the low-density limit p — > to the free-space K~N cm. forward scattering 
amplitude F K - N (p, y/s), pis the relative K~N momentum, s = (E K +E N ) 2 — 
{Pk + Pn) 2 is the Lorentz invariant Mandelstam variable equal to the square 
of the total K~N energy in the two-body cm. frame, and the nucleon 
energy is approximated by its mass in the kinematical factor in 
front of F K - N . The in-medium amplitude F K - N (p, y/s; p) in this work is a 
chirally motivated amplitude constructed within a full octet 0~ meson-octet 
l/2 + baryon coupled channel separable interaction model [13, 15] which in 



its latest next to leading order (NLO) version NLO30 [15[ incorporates the 
recent SIDDHARTA data for the atomic K~Y{ Is level shift and width j^j. In 
this separable interaction model, the in-medium coupled channel scattering 
amplitudes assume the form 

Fij(p,p' } Vs; p) = gi(p)fij(Vs; p)9j(p), (2) 
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with form factors gj(p) = a|/(p 2 + of). The momentum dependence in- 
troduced by the form factor Qk-n m the separable interaction model of 
Refs. 14|, ll5( is relatively weak for the applications discussed in the present 
work and is secondary to the strong energy dependence of the reduced am- 
plitude f k~ n generated by the A(1405) subthreshold resonance. 



2.0 

1.5 

1.0 

0.5 



NLO30 K" N average, p=0.5p 
free 

without SE 




o 50 
E-E th (MeV) 



100 



Figure 1: Energy dependence of the cm. K~N reduced amplitude §S§ in version NLO30 
of the chiral model [15[ below and above E t h — tuk + — 1432 MeV. Dashed curves: 
free-space amplitude; dot-dashed curves: Pauli blocked amplitude at 0.5po; solid curves: 
including meson and baryon self energies (SE), also at 0.5po- 



Free-space and in-medium reduced amplitudes for half nuclear matter 
density p = 0.5po are shown in Fig. [I] for the isospin-averaged combination 

1 3 1 

fK-N = -^(fK~p + fK-n) = ^//=1 + ^/ /= 0' (3) 

corresponding to symmetric nuclear matter. Similar results are obtained at 
full nuclear matter density po = 0.17 fm -3 . Of the two in- medium amplitudes 
shown in the figure, the one marked "without SE" imposes Pauli blocking on 
intermediate KN states for p ^ [8| , whereas the one marked "with SE" adds 



self consistently hadron self energies in intermediate states [16|, following a 
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procedure suggested in Ref. [17|. The real part of all three amplitudes ex- 
hibits strong energy dependence, switching from weak attraction above K~N 
threshold to strong attraction below threshold. As a rule of thumb, Re / = 1 
fm translates into a sizable attraction Re Vk- ~ —100 MeV. The imaginary 
part of these amplitudes exhibits a peak, related to the subthreshold A(1405) 
resonance, with a steep decrease at lower energies, becoming vanishingly 
small near the 7rS threshold about 100 MeV below the K~N threshold. In 
between the two limits of the energy scale in Fig. [TJ E — E t h = ±100 MeV, 
the three amplitudes differ appreciably from each other. At threshold, in 
particular, the real part of the "with SE" amplitude is about half of that 
"without SE", corresponding to a depth —Re Vk-(po) ~ 40 — 50 MeV, in 
agreement with Ramos and Oset [^j]. 

It was recognized in the early 1970s that the strong energy dependence of 
the two-body amplitude, particularly in the subthreshold region where the 
KN quasibound state A(1405) dominates, provides the underlying structure 
for K~ nuclear interactions at and near the K~ nucleus threshold 18l. ll9l 
This idea has been reformulated and applied recently in Refs. 12|, |l3 
a comprehensive study of kaonic atoms. The essential idea is to replace 
the two-body variables p and a/s of the in-medium scattering amplitude by 
appropriate density dependent averages in the nuclear medium. This may 
be summarized by the following relationships: 

-> E th - B N - B K - £ N p- - £k^, (4) 

ZrriN ZrriK 

upon neglecting quadratic terms in the binding energies B K = m K —E K , B N = 
m N — E N near threshold (E th = m N + m K ), and 

V\ P' 2 -> ^NU(2m K ^- + 2m,A), (5) 

2771 TV 2771k 

where £n(k) = rnN(K)/(?nN + ttik) i n both of these substitutions. Replacing 
in Eqs. (@| and (jSJ) the kinetic energy p k /(2ttik) in the local density ap- 
proximation by — Bk — Re Vk- (p) where Vk- = V K - + V c , with V c the K~ 
Coulomb potential generated by the finite-size nuclear charge distribution, 
and approximating the nucleon kinetic energy p 2 N /{2m^) in the Fermi gas 
model by 23 (p/p ) 2/3 MeV, Eqs. ® and © become 

v ^^ J E th - J B JV -^ J B^-15.1(-^) 2/3 + ^ReV^-(p), (in MeV) (6) 

Po 
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where all the terms following E& on the r.h.s. are negative, thus imple- 
menting the anticipated downward energy shift into the K~N subthreshold 
energy region, and 

P 2 ~ ^K[2m K 23(p/p ) 2 / 3 - 2m N (B K + Re V K -(j>))], (in MeV) (7) 

where both terms on the r.h.s. are positive for attractive potentials Vk-- The 
dominant contribution in Eqs. (jB]) and fl7]) arises from ReV^-(p), resulting in 
a downward energy shift of up to 60 MeV as demonstrated in Fig. [21 and in 
values of p(po),p'(po) as high as 275 MeV/c for K~ nuclear potential depths 
reaching 180 MeV in phenomenological studies jlj]. These momenta are well 
within the NLO30 momentum dependence scale a^ N = 700 MeV/c. 
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Figure 2: Subthreshold energies as function of nuclear density, see text. 



Having transformed the dependence of the in-medium scattering ampli- 
tudes Fx-jip, y/s, p) (j = p,n,N) on p and y/s into a density dependence, 
we denote the resultant in-medium scattering amplitudes by J-'x-jip)- In 
order to allow for different proton and neutron distributions in the actual 
calculations detailed below, the in-medium amplitude J-'k-n(p) which should 
substitute for F K - N (p, y/s, p) in the construction of the -^"-nucleus potential 
Eq. (CQ) is further replaced by an effective in-medium amplitude J-^.^p): 

J r K-N(p)p( r ) = J r K^p{p)pp{r) + J 7 K-n{p)pn{r), (8) 
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with p p and p n normalized to Z and N, respectively, and Z + N = A. The 
reduced amplitudes fx- P and fx-n are evaluated at yfs given by Eq. (jBJ), 
where the K~ atomic binding energy Bk is neglected with respect to the 
average nucleon binding energy B^ ~ 8.5 MeV a A similar approximation 
is made in Eq. (J7J) for p 2 when using the form factors gx-N^p) of Eq. (j2J). 
The i^~-nucleus potential Vjf- (p) is calculated by requiring self consistency 
in solving Eq. (jHJ) with respect to Re Vk-, i.e., the value of Re Vjc-(p) in 
the expression for yfs and in the form factors gx-N has to agree with the 
resulting Re Vk-(p)- This is done at each radial point and for every target 
nucleus in the data base. 

3. /^"-nucleus scattering amplitudes 

The present model transforms the energy dependence of subthreshold 
effective amplitudes, Eq. (jSJ), into density dependence. This transformation 
is hardly sensitive to the nucleus involved, as is seen in Fig. [2] calculated in 
the "without SE" version of the NLO30 model. In the "with SE" version 
the energies for a given density differ from the plotted values by 2-3 MeV. 
The energy shifts do not vanish for zero density because we have used a 
fixed average nucleon binding energy of 8.5 MeV. Replacing it by a position- 
dependent B N — > B N p{r)/p, in order to satisfy the low-density limit, causes 
the energy shift to vanish far outside the nucleus, with minor overall effects 
on the present results. It is seen from the figure that, e.g., for a density 
of 50% of nuclear matter density, the downward energy shift is ~ 40 MeV 
which, from Fig. [H implies a real amplitude at least twice larger than the 
threshold value in the "with SE" version. 

The density dependencies of the real and of the imaginary part of the 
effective amplitude F'k-n are of particular interest because they are related to 
characteristic features of phenomenological optical potentials. It was shown 
already in 1993 [7] that with empirical density-dependent potentials, where 
the effective K~N scattering amplitude within a tp model depends on the 
density, improved fits to the data were obtained compared to fits using fixed 
amplitudes. It was observed that in addition to the increased depth of the 



lr The precise value used for Bn in our kaonic atom global fit hardly matters within rea- 
sonable limits. Sensitivity to Bn is ex pec ted in studies limited to light nuclei, as exhibited 
recently by analyzing FINUDA data 2l| of A hypernuclear formation with stopped K~ 
mesons on targets from lithium to oxygen 22 1 . 
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best-fit real potentials, these were characterised by compression relative to 
the corresponding nuclear densities, with r.m.s. radii of the real potential 
smaller than the corresponding nuclear radii. Reduced r.m.s. radii of optical 
potentials mean that the underlying in- medium K~N interaction increases 
with density, a robust feature that is insensitive to details. 
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Figure 3: Effective amplitudes as function of nuclear density in model NLO30 without SE. 




Figure [3] shows the NLO30 effective amplitudes for Ni and Pb as func- 
tion of nuclear density, calculated in the "without SE" version. Qualitatively 



similar results are obtained also in the "with SE" version, see Ref. [12j. Re- 
gions of low density, i.e. large radii, are the most effective in determining 
the r.m.s. radius of a distribution and the sharp rise of the real part in the 
extreme surface region can lead to compression of the real potential. The 
opposite dependence is observed for the imaginary part, thus implying infla- 
tion of the imaginary potential relative to the nuclear density. Quantitatively, 
however, the change of r.m.s. radii relative to nuclear densities is found to 
disagree with the empirical trends of Ref. p}, with too little compression for 
the real part and far too strong inflation for the imaginary part. Inevitably 
this is reflected in the quality of agreement with the data, as demonstrated 
in the next section. 

Strong-interaction effects in kaonic atoms are dominated by absorption, 
as is evident from level widths being significantly larger than the correspond- 
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ing level shifts. Moreover, the shifts are always repulsive although the real 
potential is attractive, again pointing to the dominance of absorption. There- 
fore it is argued that the above marked decrease of the imaginary part of the 
effective scattering amplitude is the main deficiency of the present model, a 
decrease originating in the sharp decrease of the imaginary part of the free 
amplitude (Fig. [T]) towards the 7rS threshold, which is typical of the single- 
nucleon approach. Consequently it is reasonable to expect that additional, 
multi-nucleon terms are required to obtain good fit to the data. 

4. /^"-nucleus optical potentials 

Strong interaction level shifts and widths in kaonic atoms have been 
calculated by solving a Klein-Gordon equation [l[ with the optical poten- 
tial of Eq. ([T]) transformed to the if~-nucleus cm. system, and where the 
in- medium K~N scattering amplitude is given by the effective amplitude 
Eq. ©: 

Vk- = --(1 + ^^Ff-Mpir), (9) 
[L A in N 

with \i the kaon- nucleus reduced mass and p = p p +p n - Two-parameter Fermi 
distributions were used for both densities, with p p obtained from the known 
charge distribution by unfolding the finite size of the charge of the proton. 
For p n averages of the 'skin' and 'halo' forms of Ref. [23] were adopted with 
the difference between r.m.s. radii given by r n — r p = (N — Z)/A — 0.035 fm. 

Figures H] and [5] show optical potentials for K~-Ni and for K~-Pb, respec- 
tively. The potentials marked NLO30 follow directly from the in-medium 
K~N amplitudes in model NLO30, without any adjustable parameters. The 
agreement with the full data set of 65 points, covering the whole periodic ta- 
ble, is poor, with x 2 per point in the range of 10 — 12. This is not surprising 
in view of the obvious deficiency of the single-nucleon approach where the 
imaginary part of the amplitude goes down rapidly towards far subthreshold 
energies. Adding to the potential an empirical term linear in the nuclear 
density does not improve much the fit to kaonic atoms data and only by fur- 
ther addition of a p 2 /po or a p{p/po) 2 term good fits to the data are possible. 
Both imaginary and real parts of the additional phenomenological poten- 
tial are then found to be dominated by p 2 or p 3 terms which are likely to 
represent multi-nucleon absorptive and dispersive contributions, respectively. 
Figures S] and [5] also show potentials obtained when adjustable bp + Bp 2 /p 
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Figure 4: K nuclear potentials for K atoms of Ni. Dashed curves: derived from in- 
medium NLO30 amplitudes; solid curves: plus phenomenological terms from global fits. 



terms are added and the four parameters b and B are determined by requir- 
ing best fit to the data. The quality of the fits is then quite acceptable, with 
X 2 per point of 2 to 2.3. Qualitatively similar results to those displayed in 
Figs. S] and are obtained also in the "with SE" version, in agreement with 



the discussion for Ni in Ref. 12 



Considering values of the potentials at the nuclear center, the additional 
phenomenological part appears too large to be regarded as a correction term 
to the basic NLO30 amplitude. However, values of the potential at the cen- 
ter are rather meaningless in the context of kaonic atom observables. The 
sensitivity of calculated level shifts and widths to the K~ nuclear poten- 
tials was found 24| to be around the nuclear surface and certainly not at 
the center. With that in mind and focusing on the imaginary potential as 
noted above, it is remarkable that the additional term modifies the shape 
of the imaginary potentials in the surface region, bringing their r.m.s. radii 
closer to empirical values. Changes of the imaginary potentials near the 
surface due to the phenomenological term are of the order of 30% of the 
NLO30 potentials, consistent with the fraction of multi-nucleon absorption 
estimated from experiments with emulsions and bubble chambers 25fl- The 
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Figure 5: K nuclear potentials for K atoms of Pb. Dashed curves: derived from in- 
medium NLO30 amplitudes; solid curves: plus phenomenological terms from global fits. 



emerging phenomenology is similar to that for V^- in pionic atom studies 
where theoretically motivated single-nucleon contributions are supplemented 
by phenomenological p 2 terms representing tcNN processes [26|. (See also 
Ref. Q). 

5. Summary and conclusions 

A simple ansatz for transforming the strong energy dependence of sub- 
threshold K~N scattering amplitudes in the nuclear medium to appropriate 
density dependent averages was presented and employed in global analyses of 



kaonic atom data, following Refs. [12|, |13|. With chirally motivated coupled 



channel separable interaction scattering amplitudes in model NLO30 [15 
that respect the low energy KN data, including the recent SIDDHARTA 
results for kaonic hydrogen jBJ, the connection between this model and deep 
real optical potentials was re-established. Effective -^"-nucleus amplitudes 
were derived self-consistently and were critically reviewed with respect to 
empirical features of phenomenological optical potentials. We focused in the 
present update on the in-medium effects arising exclusively from the strong 
energy dependence of subthreshold K~N amplitudes, using for this purpose 
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the "without SE" version of the NLO30 in-medium model. The introduction 
of self-energy effects in the "with SE" version is necessarily model dependent 
to some extent. Nevertheless, all of our findings and conclusions hold true 
in both "without SE" and "with SE" versions, with minor differences exhib- 



ited already within earlier versions [13|. In the present update, as well as in 



the preceding studies [12j, |13| , the steep decrease of the imaginary part of the 
amplitude as function of the nuclear density, due to the single-nucleon nature 
of the model, was identified as a major deficiency of the single-nucleon ap- 
proach. This conclusion is valid also upon adding effective .^"-nucleus ampli- 
tudes generated by the p-wave £(1385) subthreshold resonance, as discussed 



in Ref. [13j where p-wave effects were found secondary to A(1405)-dominated 
s-wave effects. Good agreement with experiment was achieved by adding to 
the potential a phenomenological part which was found to be dominated by 
a p 2 or a p 3 term. Including systematically multi-nucleon processes should 
be the next step in trying to obtain .^"-nucleus potentials from in-medium 
K~N interaction input. 
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